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Abstract 

We calculate the finite vacuum energy density of the scalar and electromagnetic fields inside 
a Casimir apparatus made up of two conducting parallel plates in a general weak gravitational 
field. The metric of the weak gravitational field has a small deviation from flat spacetime inside 
the apparatus and we find it by expanding the metric in terms of small parameters of the weak 
background. We show that the found metric can be transformed via a gauge transformation to 
the Fermi metric. We solve the Klein-Gordon equation exactly and find mode frequencies in Fermi 
spacetime. Using the fact that the electromagnetic field can be represented by two scalar fields in 
the Fermi spacetime, we find general formulas for the energy density and mode frequencies of the 
electromagnetic field. Some well-known weak backgrounds are examined and consistency of the 
results with the literature is shown. 
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I. INTRODUCTION 


The quantum vacuum is a fundamental concept in theoretical physics and its properties 
has been widely investigated in the literature of quantum gravity and string theory. The 
theory of quantum helds in curved spacetime which is believed to be the low energy limit of 
the ultimate theory of quantum gravity, has predicted famous quantum effects in the presence 
of gravity. In general, due to the lack of the global symmetries in the spacetime manifold, 
quantum effect considerations in curved spacetime are mainly limited to the analysis of the 
local quantities such as the vacuum expectation value of the energy-momentum tensor, i.e. 
< |T^i,(a;)| >, in some point x. In fact, the most famous results of the semi-classical theory 
of gravity like the Hawking radiation and the particle production in the expanding universe 
has been achieved from the analysis of the < \T^^\ > in the related curved backgrounds. The 
most famous vacuum state effect is the Casimir effect. An important aspect in the researches 
around the Casimir effect in curved spacetime is that the characteristics of the vacuum 
state are apparently dependent to the geometry of the background spacetime. We are also 
motivated to see explicitly such a dependency in this paper. Furthermore, computation of the 
energy, i.e. < |Too| >, has also been done for large number of problems in various spacetimes 
and sometimes [1] has helped us to confirm the validity of the principle of correspondence 
in the context of the Casimir effect. Finding the total gravitational force on a set of two 
conducting Casimir plates [2] is a typical example. So, due to the importance of the stress- 
energy tensor, we will consider the < \T^^\ > for the plates in a general weak background. 

The Casimir effect arises when there is a boundary in our problem and it predicts a force 
between two uncharged conducting metals in the presence of a quantum field. The effect has 
been measured to a great accuracy |3]. We use the zero-point energy approach here although 
it is possible to find the Casimir force and energy without any reference to the zero-point 
energy |3]. We may also have the Casimir effect without having a boundary at all. In fact, 
some non-trivial topologies in curved spacetime do the same job as a boundary does[n],[S]. 
The Casimir energy in curved spacetime has been also analysed by many authors ([iiEiimn] 
and references there in). Recently, a Casimir apparatus consisting of two ideal conducting 
parallel plates in the weak field limit of the Kerr and the Horava-Lifshitz spacetimes has 
been studied in [THTO]. A purpose of this paper is to generalize the above analysis for scalar 
field doing exact solution of the Klein-Gordon equation in a general weak gravitational field. 
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Also we extend the method for the case the electromagnetic held is present inside the plates. 

The organization of the paper is as follows. In section II using the fact that the apparatus 
has composed of tiny pales, we will hnd the small deviations of the metric from hat spacetime 
inside the apparatus. In fact we will expand the metric up to hrst order in terms of the 
parameters of the general weak gravitational held. In section III the Klein-Gordon equation 
will be solved exactly inside the apparatus using the metric obtained in the previous section. 
In section IV mode frequencies inside the apparatus will be obtained under the inhuence of 
both Dirichlet and Neumann boundary conditions for the scalar held. Generalization to the 
case where the electromagnetic held is present will be done using an interesting property of 
the Fermi spacetime in section V. Gomputation of the energy momentum tensor for the scalar 
held for both Neumann and Dirichlet boundary conditions are done in section VI. Also the 
electromagnetic energy density is obtained in this section. Well known weak gravitational 
helds are examined in section VII. The electromagnetic energy density for the far held limit 
of the Kerr spacetime and the Horava-Lifschitz theory of gravity are of special interest. The 
hnal section is devoted to the conclusion. 


II. TRANSFORMATION OF THE METRIC OF A WEAK GRAVITATIONAL 
FIELD INTO THE FERMI METRIC 


In the 1+3 formalism of general relativity, stationary spacetime metric is dehned by 

ds^ = goo{dx^ — Aidx^Y — df', ( 1 ) 

where Aj = — is the so called gravitomagnetic potential and 

( 2 ) 


dP = 'jijdx^dx^ = {—gij + )dx^dx^ i,j = 1,2,3., 

goo 

In the weak held slowly rotating limit (<I> << l,n << c), the metric Q is equivalent to 


, . , 2^ 2 A.V, . , , 


2<F 

2 )PdP — (1 - :^)6ijdx^dxP 


( 3 ) 


The explicit form of a general weak gravitational held line element is as follows (see (19.13) 
in [26] ) 

2M 2M‘^ 1 1 

ds^ =[1 - — + ^ + 0{^)]dP - + 0{^)Wdt 

2M 3MA , , 

— 1(1 H-1- 5-)+7 + gravitational radiation terms that die out as U[-)Ydx dx% 

rjrt '' IT 


( 4 ) 
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FIG. 1: The Casimir apparatus far from the center of the source of a weak gravitational field. Z 
axis coincides the r axis in the equatorial plane and the plates are separated with a small coordinate 
distance I «1. The boundaries are at z=0,z=l 

Comparison between (|^ and (|^ shows that $ = — is the newtonian potential, 

= ^ijkS^^ is the gravitomagnetic potential and v* = We need to have an isotropic 

coordinate representation of the metric in the next sections and equations 0,0 and 0 
shows us the way we can construct an isotropic coordinate representation of any weak grav¬ 
itational held. 

Fig. E demonstrates the apparatus in a weak gravitational held. Two ideal conducting 
plates are in a small coordinate distance I from each other. Inside the plates, the metric has 
a small variation relative to the hat spacetime metric. To hnd this variation, we adapt a 
rectangular coordinate system having the origin at one of the plates (the one which is closer 
to the source and have coordinate distance R from it) and expand the metric 0 inside the 
apparatus in the neighborhood of the point r = R . The overall size of the apparatus is so 
small that we can assume r = R + z: 

dh (t 

gpiv{.r,0) = + h^^{R + z,0) = + hf,^{R,e) H-+ ^(ez^), (5) 

in which R » z, « 1- For the case of a static spacetime, the components of the 
metric 0 can be written in the following form: 

g^iv = 1 -|- 27 -|- 2\z -f 0 ( 7 ^^), ( 6 ) 

where 7 < 1, A < 1 are constants and use is made of d* = 7 -|-A 2 ;-l- 0 ( 7 ^) , 7 = A = 
Concerning the form of A*, the above expansion satishes 1 2=0 < 1 provided that 
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6 = 71/2 i.e. in the equatorial plane. So (|^,(|^ are also valid for the case of the far held 
limit of the Kerr spacetime and we will back to it in the examples in the section VI. 

Motivated by the above discussion, in this paper we analyse the general case of the 
spacetime of the form 

ds^ = (1 + 27o + 2\Qz)dt^ ~ (1 + 27i + 2\iz)[dx^ + dy^ + dz^], (7) 


in which 7 o, Ao, 7 i, < 1 . 

To solve the Klein-Gordon equation, it is better to recast the metric Q to the known 
Fermi metric. We use the linearized weak held regime of general relativity and change the 
variables with the aid of the following gauge transformation: 


9^iv T << 1, Cii,u 

x'^ = 

Ct = (70 + -^0(2: — 2: ))t, C = (70 + -^0(2: — 2: ))t, 
Cx = -(71 + >^iz)x, C = (71 + >^iz)x, 

Cy = -(71 + ^iz)y, C = (71 + ^iz)y, 

Cz = -llZ - XlZ^, c = llZ + 


which we have assumed h',, 
explicitly we have: 


0 to force the spatial sector of the metric (7) to be hat. More 


t = t + 7 ot, 

x' = X + {ji + Xiz)x, 

y' = y + {ii + Mz)y, 

z' = Z + JiZ + ^^iz'^- 


(9) 


the metric then takes the following form up to hrst order in the parameters 70 , Ao,7i, Ai 


ds'^ = (1 + 2Xoz)dt‘^ — dx'"^ — dy'"^ — dz"^. 


( 10 ) 


The gauge transformation, however, changes our primary problem as follows. In fact ac¬ 
cording to the last equation in ( 1 ^ the boundaries must be transformed from z = 0 and 


z = I in the spacetime (^7| to . 2 ' = 0 and z' = I + 'jil + hXiP in the spacetime (10). Another 


change that the gauge transformation brings into the problem is the rescaling of time in ([^ 
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by the factor 1 + 70 . This corresponds, in turn, to dividing the mode frequencies a; by 1 + 70 
because of the presence of the factor in our solution of the Klein-Gordon equation in 
the next section. So Time must be re-inverted after solving the Klein-Gordon equation when 
we obtain the mode frequencies. The net effect of the rescaling of time is that the final mode 
frequencies must be multiplied by the factor 1 -|- 70 . We drop the dashes ' on x\y\z from 
now on using the new boundary conditions instead. 


III. EXACT SOLUTION TO THE MASSLESS KLEIN-GORDON EQUATION IN 
THE FERMI METRIC 


The massless Klein-Gordon equation is 

=0 , 0 = (11) 

Since the spacetime is spatially flat we assume the following form for the solution 

$(a;) = (12) 

where C is normalization constant determined through the commutation relations 

(<hi(x), <hj(x)) = dijS{ki - kj). (13) 

The scalar product is dehned as 

($i,$2) = -*/' $i(x)V^<l';(x)[-0s(x)]5n'®dS. (14) 

Jy. 

in which = d^j^z and dS spans the space between the plates. Under the above assumptions 
equation ( pT| reads 

(1 + 2Az)Z"(z) + AZ'(z) + (co^ - (1 + 2Az)ki)Z(z) = 0, (15) 

where ' denotes derivation with respect to z and k]_ = k^ + k^. Another variable change 
V 1 + 2 A 2 

(1 + 2Az)V"(z) + 3\V'(z) + (co^ - (1 + 2\z)kl)V(z) = 0. (16) 


'’(0 = ASc yields 


Appearance of factor 3 in front of the second term, introduces a signihcant simplihcation 


when we change the variable to T(z) = exp(k±z)V(z). It recasts (16) into 

(1 + 2Az)T"(z) + (3A - 2A;^(1 + 2Az))T'(z) -F (cj^ - 3Ak±)T(z) = 0 (17) 
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A simple reparametrization of this last equation via u = ^(1 + 2\z) end up with 

uT"{n) + (? - u)T[n) + - ?)T(,.) = 0 (18) 

This is the known Kummer’s differential equation 


uT"{u) + {B- u)T\u) - AT{u) = 0 (19) 

o 2 q 

in which A = ^ — , B = ^. Kummer’s differential equation is not suitable for next 

considerations and we transform it via T{u) = u~^e^W{u) to another known form called 
Wittaker’s differential equation: 


in which in our case /r 
set of solutions 








)W{u) = 0 


( 20 ) 


B-l 

2 


K- = B-2A 
, K — „ 


4 k I X ' 


Equation (20) has two independent 


and their asymptotic behaviour is as follows 




r(l + 2/i)e2n M —)■ cx) 

3 

Ui M —0, 


( 21 ) 




— - K 

e 2 u 


r(2M) 

r(i+/2-K) 


M2 


M —)• CX) 
M —)■ 0 


( 22 ) 


In general, is the acceptable physical solution which is hnite at inhnity. In the problem 
under consideration we must choose a linear combination of the two, due to the fact that 
both are hnite in between the plates. The exact mode functions are: 


= ti-i|2l(w,n)in,i(“) + B(w,l;i)A4 ,(u)]e-‘"‘-“-“'-“»'' 


(23) 


The asymptotic form of (23) for small value of Aq can be written as follows (see Appendix 

A) 


Uz) = Co{u, k^)igooS{z))-lstni VSdz + (24) 

Jo 

in which S = ^ — k'j_, goQ = 1 + 2Ao2;. In the next section we use this asymptotic form 
to extract mode frequencies for both Neumann and Dirichlet boundary conditions on the 
plates. 
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IV. MODE FREQUENCIES FOR NEUMANN AND DIRICHLET BOUNDARY 
CONDITIONS FOR THE SCALAR FIELD 


We put the approximation f y/Sdz ~ y/bz + , a = —2a;^Ao ,b = uj'^ — kj_ into 


(24). From the Dirichlet boundary condition = 0) = 0 we have 0o = 0 and from 
(p^iz = Z + 71 / + = 0 we have: 




yfSdz = 


nvr. 


(25) 


After careful expansion of (|25|), the mode frequencies proved to satisfy the following relation: 

(26) 


UJ' 


'{1 — Ao[Z + 7iZ + 2-^1^^]} ~ + (^ 


UTT 




Note that the factor / + 71 Z + = Ip = y/^dz = ■\/Y+~2rj\+~2J^dz is nothing 

but the proper distance between the plates and so we have 

L 


u = a;o(l + Ao-|) 


(27) 


in which ooq = yjkp + (y^)^ , u = 0 ,l, 2 ,... denotes proper (or the corresponding flat space 
) mode frequencies in the local Lorentz frame of an observer comoving with the plates. As 
stated in the previous section, the hnal mode frequencies will be obtained by multiplication 
of the factor 1 + 70 due to rescaling of time during the gauge transformation (|^. So we have 
the hnal mode frequencies inside the Casimir apparatus for the spacetime ([^: 


u! — Ci;o(l + 7 o + -^ 0 “) 


(28) 


Mode frequencies are inhuenced only by goo component of the metric from the point of view 
of a proper observer. 


The Neumann boundary condition dz(j)\z=o = 0 imposed on (24) gives 0o = |: 

46 


_n 

— \z=0 — U 

dz 


tan{ y/Sdz + cj)o)\z=o = 


4y/S{gooS) 

iidooS) 


1 2=0 


0(Ao) 


—)■ 00 


(29) 


Another Neuman boundary condition dz(t>\z=i = 0 end up with 


dz 


cot{ / y/Sdz)\^=i = 


^VSigooS )^ 

iidooS) 


46 


\z=l 


0(Ao 


—>■ 00 


(30) 


which in turn result in (26) and (27) again. 















V. GENERALIZATION OF THE FORMALISM WHEN THE ELECTROMAG¬ 
NETIC FIELD IS PRESENT INSIDE THE PLATES 


The electromagnetic field has two physical degrees of freedom and it is known that in 
the Rindler spacetime the electromagnetic field can be represented in terms of two scalar 
fields satisfying the Klein-Gordon equation separately m The photon propagator and 
the energy-momentum tensor of the electromagnetic field in a weak gravitational of the 
Fermi spacetime has also been obtained in [T3]. However in [T3], the computations has been 
done through a lengthy and cumbersome method of the green functions. In a paper by 
the author [TH], it is shown that the energy density (i.e. the 0-0 component of the energy- 
momentum tensor) of the electromagnetic field in Fermi spacetime is exactly the some of 
the energy density of the two scalar fields mentioned in [18]. The method was done without 
any address to the Green function method frequently used in the literature. Here, we briefly 
review the relationship between the two scalar fields and the electromagnetic field in Fermi 
spacetime. 

The spin one vector field in curved spacetime in the Lorentz gauge satisfies 


= 0 , V^A^ = 0 , /i = 0,l,2,3. 


(31) 


It can be shown that in Fermi metric, the Ricci tensor satisfies = O(A^) and the second 


term of the wave equation (31) must be ignored. Furthermore, because the metric (10) is 


spatially fiat, the Lorentz gauge in (31) can be broken into two independent parts m- 

V“Aa = 0 , a = 0,3 = {t, z), 

V*Ai = 0 , i = 1,2 = {x,y). 

In which 


(32) 


A = EijV’lp , = EoiVV’- 


(33) 


and 



0 

— 1 — A 2 ; 



(34) 


We know also that both and 0 satisfy the Klein-Gordon eqaution separately (see the 
appendix in [H]). Boundary condition for the electric field on the plates is £^( 2 ; = 0) = 
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Ez{z = 0) = 0 and E_l( 2 ; = 1) = Ez{z = /) = 0 which in turn can be recast into boundary- 
conditions on 'ip and 0 according to equation (33). In [TH] it has been shown that boundary 


conditions for the electric and magnetic helds return a Dirichlet boundary condition on 0 
and a Neumann boundary condition on We have shown in previous section that both of 
this conditions will be ended up to a same frequency shift. As a result mode frequencies in 


(28) are also valid for the electromagnetic held. 


VI. THE ENERGY-MOMENTUM TENSOR 


This section has three subsections. In the hrst subsection, the electromagnetic energy- 
momentum tensor will be represented in terms of the energy momentum of the scalar helds 


mentioned in (33). In the other two sections the energy momentum tensor of the scalar and 


vector helds will be calculated. 


A. The relationship between the energy-momentnm tensor of the scalar and vector 
fields 


The vacuum expectation value of the quantum energy-momentum tensor is defined as 


< 0|T^^|0 >= 

k 


(35) 


The classical energy-momentum tensors for the scalar and vector fields are 


T. 


Scalar 


= d^(t)d^(t) - dxcpdgct), 


rnvector _ rpGhost _i_ rpGauge i rpMaxwell 

fliy ^ ^ /ii/ ’ 


pGhost 

- fii/ 


= c'[A„A';^, - + ;5(3iyni. 


rjiGauge _ _ * — r* n — 

■J- ,i7j — G-z/Uriy 


A6>. 


• /ii/ 




T. 


Maxwell 

fiu 


— z9tiuE^^Exe — E^bE^ 


Ou- 


(36) 


As is well-known, in the quantum level, the contributions of ghost and gauge fields in the 
electromagnetic energy-momentum tensor cancel each other [13] and we only concern the 
Maxwell sector of the energy-momentum tensor. Expansion of in terms of scalar 
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fields V’ and 0 shows that (see Appendix B) 


rj.Maxv,ell ^ 1 ^ ^ 2Xoz)B^), 

T^axwell ^ ^ ( 37 ) 

T^^axu,eii ^ 1 _ 2\^z){E‘^ + (1 + 2\oz)B^)g,j - (1 - 2Xoz)E,E, - B,Bj, 


in which E'^ = gijE^E^, B^ = gijB"B\ gij = -6ij, Ei = Fq*, Eij = SijkB^ . Using 


E^y = Afj,^y - Ay^^ and (32) we have: 


E-rx = —f/Ca;\/l + 2Xzdz'llj + ukycf), 

Exy = —ikyy/l + 2Xzdz'4> ~ ^kx(j), 

Et-z = —Vl + 2Xzk^i/j, 

Eyz = ikxdz4> + Ci;(l — 2Xz)Vl + 2Xzky'il), 

Exz = —ikydz(j) + C(;(l + 2Xz)Vl + 2Xzkxip, 
^xy kj_(j), 

^ = {Er X, Er y, Er z), 

( ^yzi Fxzi ^xy') 


(38) 


in which we have used the general form of the wave function (12). Quadratic products of 
helds , F^Fj and BiBj produce terms like 'ipdz'i)*-, 4>4’* which have no 

contribution when the expectation value is taken because of the fact that '0 and 0 are not 
correlated and belongs to independent Hilbert spaces. We calculate 0 — 0 component of the 
energy-momentum tensor for the electromagnetic held: 


<0|r„t,|0> ^ZJi^kzli(u,^ + {l + 2\z)klM-^ + l{l + 2\z)\dA?} (39) 


< 0|r„"“ |0 > = - X) /> +(1 + 2Az) < 0|B^|0 >)} 

LO ^ 

I Ol^’i + Pr, + > +(1 + 2Az) < 0|F7 + Fl + 4,10 >)} 

= [ (fk±{k]_[-{{u'^ + {1 + 2Xz)k]_)\(()\‘^ + {1 + 2Xz)\dz(f>\‘^} ^ 

LJ ^ 

+ A A 2A2;)fc^)I'^l^ + (1 + 2Xz)\dz'ilj\^'\]] 

= fcl[<0|Tot,|0> + <0|Tot|0 >] 
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After a lengthy but straightforward calculation, other components of the energy-momentum 
tensor in both sides are related to each other as follows: 


< > = kl{< 0|Tot,|0 > + < 0|Tot|0 >}, 

< > = -kl{< 0|Tfi|0 > + < 0|Tit|0 >} + 2(1 - 2\z)kl{< 0|Tot,|0 > -f- < 0|To^o|0 >}, 

< > = -kl{< OlT^tlO > + < 0|T2t|0 >} + 2(1 - 2Xz)kl{< 0|Tot,|0 > + < 0|Tot,|0 >}, 

< 0|T3““^“’"“|0 > = kl{< OlTg'^glO > + < OlTg^glO >}, 

< 0|ToY“^“’"“|0 > = kl{< 0|To^i|0 > + < 0|Tot|0 >}, 

< > = kl{< OlTo'^alO > + < 0|Tot|0 >}, 

< > = -kl{< 0|To^3|0 > + < OlTg^glO >}, 

< > = -kl{< OlTfalO > + < 0\Tf^\0 >} + 2(1 - 2Xz){< 0 \T^q\0 > + < 0|Tot,|0 >}, 

< 0|T2^“^“’"“|0 > = -kl{< 0\Ti^\0 > + < 0\T^3\0 >}, 

< > = -kl{< OlTfglO > + < OlTi’^glO >}. 

Note that can be absorbed into C. 


B. Energy density for Dirichlet and Neumann scalar fields 


This section is devoted to the calculation of the energy density for the Casimir apparatus 
via the direct method without any reference to the traditional Green function method. Using 
the approximations 

(^oo5)"3 ~ 1 - (^ + ^)z, 


/ \/Sdz zz \/bz-\ — %=z‘^ ,a = —2u'^X ,6 = a;^ — 

J AVb 

sin{\/hz H- = sin{\/hz) H- ^z^cos{\/hz)^ 

4V& 4V6 


cos{'/bz H- 1=^"^) = cos{'/bz) — _ 

4v6 4v6 


z sin 


in{\/hz). 


We expand the wave function (24) and hnd up to hrst order in A: 

Z{z) = Zq{[1 — (— + ^)z]sin{Vbz) H— ^z‘^cos{Vbz)} .Dirichlet 
2 4o 4 a/6 

Z{z) = Zo{[l ~ (2 ~ sin{Vbz)} .Newmann 


(42) 


( 43 ) 


Zq can be absorbed also in Cq. 
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The energy density is defined by £ = < 0|T^j,|0 > where is the lapse vector 

normal to the hypersurface 2 ; = constant i.e. n^ = d^z = ^+ 2 ^ 0 ^ ^^’ mean energy 

density so has the following form: 


" = I 


< 0|T^J0 > dz 


dz 


— J d?k^- J {-(w^ + (1 + 2Xoz)k\)Z{zY + -(1 + 2\Qz)\dzZ{z)\^] ^ 2Xqz) 
= J2[ d^k±H{u, k±) 

LJ ^ 

Calculating the factor H{u,k±) for both wave functions in (45) results in 

HDirichlet^^^ =C'2{ia;2^i + ]^klF2 + ^Fs}, 

z sin{2\/bz) a .z'^ zsin{2y/bz) cos{2^/bz), 


a , z 


1 o, cos(2 


^[—sin(2Vbz) + (——z‘^) _ 

AVb^2b ^ ’ ^2b ’ 2Vb 


-]}i; 


z sin{2y/bz) a .z'^ zsin{2y/bz) cos{2y/bz)- 

- iVf, ■ * idn - 8^- 


++(4 - 


{2\/bz)- 


AVb 2b^ 


-2b 


2Vb 


(45) 


.z sin(2\/bz) a .^rZ^ zsin{2\/bz) cos{2y/bz) 

- + (1 - + l^^)cos(2Vimi 


4 '2b 


•uraan 


'{oj, / c _ l ) + 2 

z sin{2Vbz) a z"^ 

^^=^2^^7b -+ 


a .z'^ zsin{2\/bz) cos{2\/bz), 

^ ^ ^ 47 ^ ^ 86 ^ 

4^b^2b^ 276 

z sin{2\/bz) a zsin{2\/bz) cos{2^/bz) 

=b + ^ 7 ^ - +a'T + 4 V 5 ^ 

_ ^|£«„(2v^z) + (i - .2 y»i'(2767 n|f, 

^ r,r^ sin{2ybz) ,a ,xr£^ zsin{2ybz) cos{2ybz). 

^==<''2 — i 7 ^ + <^-")l 4 -^ 7 I - 


^|^sm(2£6z) + ( 


( 46 ) 
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Simplification of the terms like cos{2\/bz) and sin{2y/bz) is possible using the fact that 


in boundaries the wave function (24) must vanish and: 


/ sfSdz = nvr — $0 

Jo 


Vbz H— %=z^ = mr — $ 0 , 


(47) 


in which 2 ; takes one of two boundary values z = 0 and z = I + ■yil + Equation (47) 

result in 

a 


sin{2Vbz) = ~ *^0 ’ cos{2Vbz) ~ 1, 


and the hnal result is as follows 


H 


Neuman 


>0; 




(48) 


(49) 


2 ' 46' 

Exactly the same H is obtained for Dirichlet boundary condition. The constant C was 
defined in ([^ and can be determined simply as: 




= N { I {[1 — (2A + ^)z]sin‘^{Vbz) H— 'J^z‘^sin{2Vbz)}dz} ^ ^Dirichlet 

' r‘ a 

{ / {[1 - (2A + ^)z]cos^(Vbz) - -i=z^sin{2\/hz)}dz} ^ , Neumann 


a 


2(27r)2a; Vo 
1 


2(27r)2u; Vo 


Ax/b 


Just like H, the constant C has also the same form for both boundary conditions up to hrst 
order in A 


= 


1 I a^(V-Q(2)^_i 


(51) 


The final result for the energy-density after all is: 

u 

e = 


[ dJki 

UJ ^ 

= 2(2.)^ 

= (1 + 7o + Ao— j d^ki 


2i2TiY 

(:uo(l -F7 o + Ao^) 


(52) 


(Uo 


2(27r)2' 


in which £0 = I d‘^k± 2^1)'^ corresponding flat spacetime Casimir energy density 

^0 = [22]. 

Up to now, we have shown that for Neumann and Dirichlet boundary conditions there 
exists the relation < 0|T^p,g^^j|0 >= (l+yo + Ao^) < ^I'^oo Fiatl^ > between flat and curved 
energy density contents of the Casimir apparatus in the weak spacetime of the metric ([^. 
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C. Energy density for the electromagnetic field 


Using the first equation of (41) it is evident that the same shift in the energy density 


obtained in (52) holds also for the electromagnetic field : 


< > =< 0|T(fo|0 > + < 0|To^o|0 > 


- {< 0|^ot),/«at|0 > + < 0|^ot),/iat|0 >}(1 + 70 + -^ 0 ^) 


= {- 


TT 


h 


1440^“ + + 

h 


(53) 


TT Lp 

In the next section we will analyse some well-known weak gravitational fields and find the 
parameters 70 , 71 , Aq, Ai in each case. 


D. Notes on the divergences 

The problem of the divergences near a perfect generic conductor first studied systemat¬ 
ically by Deutsch and Candelas [2H]- They found that the energy-momentum tensor near 
the surface behaves like: 

< 0|Too|0 >= ^ + •••) (54) 

where ci, C 2 are constants and e is the distance from the surface of the ideal boundary. In case 
where there is a conformal invariance in the action ci = 0. The divergences originates from 
the unphysical nature of classical ideal conductor boundary conditions. It has been shown 
[28] that we can remove the infinities of the total energy (and not the energy-momentum 
tensor) of the plates using the zeta function regularization unless for the case that the zeta 
function has poles itself. On the other hand, the cut off regularization method suggests the 
removal of the divergences via ad hoc although in this method it still remains a logarithmic 
ambiguity [22] in the energy density. For the imperfect conductors (more realistic boundary 
conditions) we can easily remove the divergences introducing some suitable cut off frequencies 
although the boundary effect may become quit large (but finite) [2l 128]. 

Now the question is that what happens when we go to curved spacetime?. Does the surface 
divergences of the energy-momentum tensor are ignorable in the semi-classical Einstein’s 
equations?. The answer is negative according to [28]. In [29| however, the authors has found 
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a way to get ride of the surface divergences of the Einstein’s held equations for the case the 
boundary is a parallelepiped. They have used a suitable cut off along with the so-called 
Estrada-Kanwal distribution theory of asymptotics to regularize/renormalize the inhnities 
and show that the energy-momentum tensor near a plane boundary, as a source, converges 
to a consistent theory when the cutoff is removed. Remarkably, the process of curing the 
divergences in curved spacetime has been continued by Milton et al [2] where they have 
shown that most of the energy between the plates is restored near the surfaces and a part 
of it resides exactly on the plates. They hnally have shown that this energies responds to 
gravity just like any other hnite energy following the newtonian relation F = ma = —Mg. 
This is expectable as this large energies are simply a part of the total energy of the system. 

We show here that, in our case, the divergences do not present in the hrst order of 
approximation that we have used here and they appear only in the higher orders of ap¬ 
proximation. To do so, we write down the explicit structure of possible divergences in the 
Einstein’s held equations along the lines depicted in [22] • The one-loop effective action W 
for the semi-classical theory of gravity is: 

r. 1 r. A SttG , 

Rfiu T ^ 

w = j ^/^LeffXx)(Fx, (55) 

1 

Leff.{x) = - Im / drri^Gp-^Xx.x), 

2 x^x Jm'^ 

in which G^'^'{x,x) is the DeWitt-Schwinger-Feynmann’s propagator. Using the DeWitt- 
Schwinger representation of the action, the asymptotic expansion of Tg//. is as follows: 

- CXD 

Leff.{x) = -(47r)-t(—- ^), (56) 

^ i=o 


in which /i is a length scale to hx the dimensional issues. The potentially divergent part of 
the effective action (the hrst three terms) [22] reads : 

1.0 1 


Ldiv = -( 47 r) 2 {--- + -[7 + In —]}[ 


"n — 4 


g? n{n — 2) 


4m Oo — 


(n-2) 


2 m^ai -|- 02] 


(57) 


in which m is the mass of the scalar held(m=0 in our case). The coefficients oq, ai, 02 are: 

ao(a;) = 1, aAx) = -R, a2{x) = —- -DR + —R^ (58) 

J , J g , apid ^gg ap g 72 
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in which Ra/ 3 -fS is the Riemann’s curvature tensor. The total gravitational lagrangian density 
can be shown to have the form (see (6.49) in [22]): 

1 


=5^ = -(R + 




SttG 


B 


A = 


B = 


IQttGb 

Am'^ 




1 1 Tn‘^ 

(4^)*n(n-2)*:rVl + 2l'' + ‘'V'’ 


m 


1 1 TTl^ 


,2 


(59) 


However, in the massless case of ours, the only non-vanishing potentially ultraviolent term 


is the one related to 02 in (52) and A, B in (52) vanish (see (6.101) in [22|)- Our calculation 
for the metric ([^ shows that 02 is of second order of approximation: 

AqAi 


R 1212 — -R 1313 — 

Ri414 = 


4(1 + 271 + 2Aiz) 


X 2 


+ 


= 0 (A^), 

AqAi 


4(1 -h 270 -h 2Ao^) 4(1 271 2Ai;2) 

A? 


— 2i?2323 — -R 2424 — R.3434 — 


L 3434 


= 0(Ag), 
= 0(A?), 


(60) 


2(1 -|- 271 -|- 2Aiz) 

Rn = R 22 = R 33 - 0{Xl), R 44 ^ 0(A2), R ~ 0(A2), DR ~ 0(A2). 


So the total lagrangian density (59) reduces to the the standard bare density 


WttGb 


R. In 


conclusion, the potentially divergent term 02 vanishes within the hrst order of approximation. 


VII. EXAMPLES:FINDING COEFFICIENTS 70,71, Ao,Ai 

In this section number of spacetimes are investigated and the parameters appeared in 
,( [53| ) will be found. To this end, we will try to hnd their weak held form according 
to and ([^. 

A. Electromagnetic Casimir energy density for the far field limit of the Kerr 
spacetime 

Recently Bezerra et al [7] studied the Apparatus for the scalar helds in the weak held 
limit of the Kerr spacetime in the equatorial plane. In that work the apparatus co-rotates 
with the local angular velocity of the spacetime i.e. the measurements had been assumed 
to be done in the point of view of a zero angular momentum observer(ZAMO). They found 


(28),(52 
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the metric inside the apparatus through some two stage successive approximation method 
as follows: 

ds^ + 2b^o)dt‘^ ~ (1 ~ 2<ho)[(ia;^ + dy"^ + dz'^] 0 — th order of approximation, 

(61) 

ds^ ~(1 + 2b\z)dt^ ~ (1 ~ 2\z)[dx‘^ + dy'^ + dz'^] first order of approximation. 


in which 6 = 1 — 2aflQ, <ho = A = — cl is the angular momentum per mass 

and flo the local angular velocity of the Kerr spacetime. The authors set <l>o = 0 for the hrst 
order of approximations. This does not works because A is related to $0 through 7 = —i?$o- 
Putting one of them equal to zero forces the other one to vanish also. Evidently, in the hrst 
order of approximations we must keep both $0 and 7 and the metric must be written as 


follows instead of (61): 


ds‘^ zedf^ — [dx^ + dy^ + dz'^] 0 — th order of approximation, 


(62) 


ds^ ~[1 + 2fe(<ho + '')z)]dn — [1 — 2(<l)o + 72 ;)]+ dy^ + dz^] first order of approximation. 


Taking the above comment and the metric (62) into account, we hnd 70 = b^o, Aq = 
67 = — 7 i = *^* 0 ; Ai = 7 = —bR^Q and arrive at the following relations for the 
electromagnetic Casimir energy density and the mode frequencies: 


'yl p 

u = a;o[l + 6<ho + 


< > = 




(63) 


,. [1 + 6<ho + b- —1- 
720/4 L ^ ow 2 J 


The scalar held Casimir energy density inside the apparatus which has been sketched in 
eq.(32) in [7] must also be corrected as follows: 


U} = + bia + 

< oir‘7-10 > = - + b^]- 


(64) 


Far field limit of the Schwarzschild spacetime is also covered by setting 6 = 1 


B. The Fermi spacetime 

The Fermi spacetime is described by the following metric 

ds^ = (1 + 2az)dt^ — dx^ — dy'^ — dz'^ (65) 
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The importance of this metric is that it is traditionally recognized as the spacetime of a 
static accelerating observer near the surface of the source of a constant gravitational held 
j. Comparison of this metric to the general metric in Q gives 70 = 0, Aq = a, 71 = Ai = 0: 


_ n I (dp, 
(U — (Uo[l + 


< > = 


.—[1 + ^1 
720/4^ 2 


( 66 ) 


< 0 | To " o “^“"|0 > = 


TT 


-[ 1 + — 1 . 
1440/4 ^ 2 ^ 


The above energy densities are exactly the results (5.2) in [I3], (3.4) in [16] and (5.4) in [T5] . 


C. The Hofava-Lifshitz gravity with a cosmological constant 


The Hofava-Lifshitz (HL) gravity is a renormalizable theory of gravity that is invariant 
under the Lifshitz scaling transformation x —&x, t —)■ Vt. This transformations manifestly 
break the space and time covariance. The anisotropy between space and time, in turn, may 
affects the Casimir effect as well. It is interesting to investigate the vacuum characteristics of 
the theory. Recently, the effect of the HL theory on the Casimir energy of the apparatus has 
been studied in [25] . The authors recommended to set a constraint on spacetime anisotropies 
in such a way that the Casimir energy modifications remain within the experimental bounds. 
Recently in 0 . CD] the same problem considered in curved spacetime in the context of 
a spherical symmetric solution of the HL theory. Finite temperature Casimir energy in 
spacetime ([^ has been analysed by the author in [36] and the following weak field limit for 
the HL theory has been calculated: 


M 3M^ 


ds — {1 -f [ 2 — -f 


M 2 


R 2R2 2cuR4 


^ M 3M2 M\ ^ . 


R 2 


, , M m2 


\— 


2R^ 2u^- 

M 2 m 2 


(67) 


]z}{dx‘^ -f dy‘^ + dz‘^) 


4R3 4wR5J 

where M = M(l-|-A). This spacetime is the weak held limit of the Park’s spherical symmetric 
solution to the IR limit of the HL theory in the presence of a cosmological constant ly [TO] . 
A may have the same role as the cosmological constant but not necessarily being a small 
parameter and a; is a constant frequently used to regulate the UV limit of the HL theory. 
Comparison between (67) and the general metric ([^ shows: 


7o 


M(l + ^) 

R 


! Aq 


M(l + ^) M(l + ^) 


R 2 


2R 


2 R 2 


( 68 ) 
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Another spacetime, which is a solution to the HL theory without cosmological constant is 
the Kehagias-Sfetsos (KS) solution. This spacetime has been discussed in [U] to obtain the 
energy density of the apparatus and is as follows: 


= fxsdt^ — fxsdp^ ~ 


(69) 


where fxs = 1 + — \/(l + ^))- P is a radial coordinate and u is the free parameter 

of the HL theory. Putting A = 0 in the Park’s solution recover the KS solution and so the 


coefficients in equation (68) are also valid for the KS solution. 


VIII. CONCLUSION 

We analysed the energy density of a Casimir apparatus consisting of two nearby con¬ 
ducting parallel plates in a general weak gravitational held. The metric in the equation ([^ 
denotes the deviation of the weak gravitational held from hat spacetime inside the appara¬ 
tus. We transformed the metric ([^ through a gauge transformation into the Fermi metric 
and then solved the Klein-Gordon equation exactly. The mode frequencies were found for 
the scalar held inside the apparatus for both Neumann and Dirichlet boundary conditions in 
terms of the weak gravitational held parameters 7 o, 7 i,Ao,Ai. This result was shown to be 
valid also for the electromagnetic held in section V. The energy density of the apparatus was 
found for both scalar and electromagnetic helds in terms of the weak held parameters. Some 
examples of weak gravitational helds were analysed in section VII. Specially the electromag¬ 
netic energy density and mode frequencies in the far held limit of the Kerr spacetime in its 
equatorial plane were obtained. The weak held limit of the Hofava-Lifshitz gravity with a 
cosmological constant was also investigated and the weak held parameters were sketched. 
Consistency of the results with the literature was checked by considering the Fermi metric. 
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Appendix A: Asymptotic form of the wave function 


In this section we find an explicit and simple asymptotic form for the wave function. As is 
apparent we need to have an asymptotic expansion with both argument and first parameter 
being large. Whittaker functions have such an expansion in terms of Airy functions [M] : 


OO 


n=0 ^ ^ 


A„(C) , A^'[(4«:)iC]>^5n(C) 


+ 


( 4 k 


n=0 ^ ’ 


}(A1) 


where C is defined as 


-(—C)2=cos ^(s/x) — Vx — x‘^ (A2) 

3 

and in our case 4kx = ^goo and so 

^2 u2 2 

X = -^goo = ^(1 + 2Az) , X < 1 or S = — -(A3) 

UJ^ goo 

If K —)■ OO then both the argument and the first parameter go to infinity. As k oc the 

4 

second term in the bracket is of order As and must be ignored to stay within the first order 
of expansion in A. Also the summation in the first term reduces to only n = 0 term and 
Ao(C) = constant. 

X -L 

The argument of the above Airy function can be written as 


(4 k)3 ^ = —{3fi:[cos ^(y/x) — Vx — x^]} 3 = —v(x) 


u —)■ OO 


(AS) 


The Airy function with large argument is as follows (section 9.7 from [30] ) 


Ai(-v)^v ^cos(|ui - ^)^(-l) 


0-2n 


. / 2 3 TT ^ 

+ sin{-v3 —- 


-3- 4-^' • (4„S)2»+i 




Ci2n+1 


}(A6) 


n=0 VS ' n=0 V 3' 

As = 0(A“^), the sm()-part and n > 0 in the cos()-part must be ignored also. The 


observation that 

cos~^(y/x) — yj X — x'^ = 


1 — X 


X 


dx = —-— / y/Sdx = — — / y/Sdz, (A7) 


2k 


and 


/ lo-i 

(-) I OC 4 

X — 1 


(AS) 
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moves the situation forward significantly: 


fK,/. oc S 4cos( / y/Sdz + 0o) , 4>o = 4 >o{k,u) 


(A9) 


A same process can be applied to the other Whittaker’s function except that according 


to [35] instead of the a argument (4k) 3 ^ in (A5) we have (4k) 3ea;p(±^)C. This difference 
changes the above process as follows: 


C—>-e*''3'(^ ^ e±"3A; , ^ =—^v 


(AlO) 


OC S' 4cos( / \^dz + (j)i) , (/>! = 0i(k, cj) 


As a result the total wave function in (22) when A —)■ 0 reads 


Mu) = Co{u, k^)igooS{z))-htn{ / ^/Sdz + 


(All) 


(A12) 


where 02 = 02(i^, A, 5, 0o, 0i) and A,B came from (22) and the following relation has 
been used: 


Acos( / \/~Sdz + 0o) + Bcos{ / \/~Sdz + 0i) = Csin{ / \/~Sdz + ^ 2 ) 


(A13) 


Appendix B: Computation of the components of the energy-momentum tensor 


In this section we hnd T*h Other components will be hnd by similar techniques. We 
do the computations in the Fermi spacetime i.e. 5^00 = 1 + = 0,5'ij = —dij. Greece 


indices run from 0 to 3 and Latins from 1 to 3. From (36) we have 


■^Maxwell ~ ^ ^ ^6 


e- 


(Bl) 


The hrst term has the following form: 

F^^F,^ = g^^g'^^F^,e,,kB^ = ... = -2B\ 


(B2) 
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where we have used SijkSiji = 26ki- The second term in (Bl) simplihes as follows: 


= gooE\ FI = 

= g^g^^Fik = g^^g^hlknB^, F^ = g^^F^m = g^^eambB\ 

F^^F^^ = -g^^g^^emabSminB^B^ 

= -fg^^[5alhn - 5aJu]B'^B^ 

= ... = -g^^B^ + B^B^. 


Based on (B1)-(B3) we hnd 


1 


T%a.u,eu = + B^) - g^^E^E^ - B^BK 


After lowering the indices we have hnally: 


T^axweii ^ _g..(^gOOF^ + B^) - g^^EiEj - BiBp 


‘O 


Now we hnd Tq*. From (36) we have 


T| 


00 


— -flUUpAyp _ pwpv 
Maxwell ~ 4" ^ ^ -^6’ 


(B3) 


(B4) 


(B6) 


(B6) 


in which the hrst term already has been obtained in (B2). The second term has been also 


obtained in (B3) and after lowering the indices again we have 

1 


E 


Maxwell 

00 




(B7) 


Now we hnd Tq*. From (36) we have 
1 


T 


•Oi _ * w pAy p puy 

Maxwell ^ -^A0 n r g 


= -F^^F) = -gooE^g^^Fmj = -goog^^SmjkE^B^ = -go^g^\E. x B)™, 
from which, after lowering the indices again, we hnd 


(B8) 


T Maxwell 
n 


■Oi 


= -(E X B), 


(B9) 
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